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ABSTRACT 

We report the results of a comprehensive study of the relationship between galaxy size, 
stellar mass and specific star-formation rate (sSFR) at rcdshifts 1.3 < z < 1.5. Based 
on a mass complete (M* > 6x 10 10 M©), spectroscopic sample from the UKIDSS Ultra- 
deep Survey (UDS), with accurate stellar- mass measurements derived from spectro- 
photometric fitting, we find that at z ~ 1.4 the location of massive galaxies on the size- 
mass plane is determined primarily by their sSFR. At this epoch we find that massive 
galaxies which are passive (sSFR ^ 0.1 Gyr -1 ) follow a tight size-mass relation, with 
half-light radii a factor f g = 2.4 ±0.2 smaller than their local counterparts. Moreover, 
amongst the passive sub-sample we find no evidence that the off-set from the local size- 
mass relation is a function of stellar population age. In contrast, we find that massive 
star-forming galaxies at this epoch lie closer to the local late-type size-mass relation 
and are only a factor f g — 1.6 ± 0.2 smaller than observed locally. Based on a sub- 
sample with dynamical mass estimates, which consists of both passive and star-forming 
objects, we also derive an independent estimate of f g = 2.3±0.3 for the typical growth 
in half-light radius between z ~ 1.4 and the present day. Focusing on the passive sub- 
sample, we conclude that to produce the necessary evolution predominantly via major 
mergers would require an unfeasible number of merger events and over populate the 
high-mass end of the local stellar mass function. In contrast, we find that a scenario 
in which mass accretion is dominated by minor mergers can comfortably produce 
the necessary evolution, whereby an increase in stellar mass of only a factor of ~ 2, 
accompanied by size growth of a factor of ~ 3.5, is required to reconcile the size- 
mass relation at z ~ 1.4 with that observed locally. Finally, we note that a significant 
fraction (44% ±12%) of the passive galaxies in our sample have a disk-like morphology, 
providing additional evidence that separate physical processes are responsible for the 
quenching of star-formation and morphological transformation in massive galaxies. 
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1 INTRODUCTION 

Following the introduction of near-IR selected galaxy sur- 
veys nearly a decade ago, it became clear that a substan- 
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tial population of massive, apparently evolved, galaxies was 
already in place by z ~ 1 (e.g. Cimatti et al. 2002; Glaze- 
brook ct al. 2004), a result which was completely unexpected 
within the context of contemporary models of galaxy for- 
mation and evolution (e.g. Kauffmann & Haehnelt 2000). 
Moreover, in line with our current understanding of galaxy 
"downsizing" , it was quickly realised that the high-mass end 
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of the galaxy stellar mass function was assembled much ear- 
lier than had been expected, with studies of near-IR selected 
samples confirming that ^ 75% of the local number density 
of M* ^ 3 x 10 11 M galaxies was already in place by z ~ 1 
(e.g. Caputi et al. 2005). 

Although the latest generation of galaxy evolution mod- 
els are now more successful at reproducing the observed 
number densities of high-mass galaxies as a function of red- 
shift (e.g. Cirasuolo et al. 2010), it has been known for sev- 
eral years that the physical properties of high-mass galaxies 
at z ^ 1 are substantially different from their low-redshift 
counterparts. Specifically, based on HST imaging of near-IR 
selected galaxies in the Hubble Ultra Deep Field, Daddi et 
al. (2005) noted that many high mass (M^ > 10 11 Mq) early- 
type galaxies (ETGs) at z ^ 1.5 displayed substantially 
smaller half-light radii than equivalently massive galaxies 
observed in the local Universe. Over the last few years, the 
size evolution of massive ETGs has been extensively dis- 
cussed in the literature and many studies (e.g. Trujillo et al. 
2007; Cimatti et al. 2008; Saracco et al. 2009; Williams et 
al. 2010) have reported that massive red/passive ETGs at 
1 < z < 3 are a factor of 2 — 5 smaller than ETGs of the 
same mass in the SDSS (Shen et al. 2003). 

Despite initial concerns that results regarding the com- 
pactness of high-redshift passive galaxies could be affected 
by surface-brightness dimming, morphological k-corrections 
or even unresolved AGN components (e.g. Daddi et al. 
2005), deep, high-resolution, near-IR imaging has con- 
clusively demonstrated the compactness of many massive 
(M* ^ 10 11 Mq) galaxies at z ^ 1 (e.g. Buitrago et al. 
2008; Damjanov et al. 2009). In addition, it is now well es- 
tablished that, within a given redshift and stellar mass in- 
terval, early- type/passive galaxies are more compact than 
late- type/star-forming galaxies (e.g. Zirm et al. 2007; Toft 
et al. 2009; Williams et al. 2010; Wuyts et al. 2011). 

In reality, the full picture is undoubtedly more com- 
plicated still, with some studies (e.g. Mancini et al. 2010) 
identifying examples of massive z ^ 1 ETGs with half-light 
radii fully consistent with the local size-mass relation and 
others (e.g. Saracco et al. 2009) finding that the popula- 
tion of high-redshift ETGs is comprised of a mixture of nor- 
mal and compact galaxies. Moreover, it is still not currently 
clear whether the compact nature of some massive ETGs at 
high redshift is related to the age of their stellar population 
(Saracco et al. 2009), or not (Trujillo, Ferreras & de La Rosa 
2011). Finally, there have been concerns raised in the liter- 
ature that comparing massive high-redshift ETGs with the 
local galaxy population is subject to severe progenitor bias. 
For example, it has been claimed that the number densities 
of compact ETGs at high redshift are actually compatible 
with compact cluster galaxies seen locally (Valentinuzzi et 
al. 2010), although this has been disputed by Taylor et al. 
(2010). Significantly, van Dokkum et al. (2010) attempted to 
circumvent the progenitor bias problem by studying the size 
evolution of galaxies selected to have a constant number den- 
sity of n ~ 2 x 10~ 4 Mpc~ 3 (equivalent to M* ~ 3 x 10 11 M 
at z = 0), finding strong evidence for size evolution of the 
form R e oc (1 + z)~ l z . 

It is thus clear that some physical process (or processes) 
must be in place at high redshift which is capable of sub- 
stantially increasing the half-light radii of a significant sub- 
set of massive high-redshift galaxies, without requiring an 



unphysical increase in their stellar masses. A widely advo- 
cated solution to the problem has been the "dry-merger" sce- 
nario, whereby dense high-redshift galaxies are transformed 
via a series of dissipationless (i.e. dry) mergers (e.g. Cimatti, 
Nipoti & Cassata 2012). Such dry mergers are typically en- 
visioned to be major-mergers (i.e. mass ratio ^ 1:3) in which 
the galaxy half-light radius increases in direct proportion to 
the accreted stellar mass. Although this scenario has several 
desirable characteristics, it also suffers from several obvious 
problems (see Nipoti et al. 2009; 2012). Principal among 
these is the prediction from N-body simulations that mas- 
sive M* ^ 10 11 Mq galaxies at z ~ 1.5 are only likely to 
undergo one major merger over the proceeding ~ 9 Gyr 
(e.g. Hopkins et al. 2010), strongly suggesting that major 
mergers alone cannot explain the size evolution of massive 
high-redshift galaxies. 

An alternative mechanism invokes size growth via a 
sequence of minor mergers, with typical mass ratios of 
1:10. The advantage of this mechanism is that in a mi- 
nor merger, radius growth is expected to be proportional 
to the square of the accreted mass (e.g. Naab et al. 2009; 
Bluck et al. 2012), allowing rapid size growth without ex- 
cessive growth in stellar mass. The minor-merger scenario 
has recently gained support from hydro-dynamical simula- 
tions (e.g. Johansson, Naab & Ostriker 2012; Oser et al. 
2012) and the observational results from van Dokkum et al. 
(2010) which demonstrate that, at a constant number den- 
sity of n ~ 2 x 10~ 4 Mpc~ 3 , the size evolution of galaxies 
is consistent with an "inside-out" growth scenario. Finally, 
it has also been suggested (e.g. Fan et al. 2008) that feed- 
back from active galactic nuclei could expel material from 
the central regions of high-redshift galaxies, thereby allow- 
ing them to expand significantly at a constant stellar mass. 
However, the increase in stellar-velocity dispersion with red- 
shift predicted by this model may be in conflict with the 
latest available data (e.g. Trujillo et al. 2011). 

In this paper we analyse the sizes, stellar masses, mor- 
phologies and specific star-formation rates of a unique sam- 
ple of massive galaxies from the UKIDSS Ultra-Deep Sur- 
vey (UDS) with spectroscopic redshifts in the range 1.3 < 
z < 1.5. Based on this sample we investigate the evolu- 
tion in the galaxy size-mass relation over a look-back time 
of 9 Gyr, using accurate stellar-mass measurements derived 
from full spectro-photometric fitting of the available spec- 
tra and multi-wavelength photometry. Armed with accurate 
stellar-mass measurements and star-formation rates based 
on a range of empirical indicators, we investigate whether 
any off-set from the local size-mass scaling relations is pri- 
marily driven by morphology, stellar population age or spe- 
cific star-formation rate. 

The structure of this paper is as follows. In Section 
2 we describe the relevant spectroscopic and photometric 
data, before describing the methods adopted to estimate 
the stellar masses and galaxy radii in Section 3. In Section 
4 we present our main results on the stellar mass - size re- 
lation and investigate the dynamical mass - size relation for 
a small sub-sample of objects with reliable stellar-velocity 
dispersion measurements. In Section 5 we outline a plausi- 
ble evolutionary scenario which is capable of reconciling the 
size-mass relation of passive galaxies at z ~ 1.4 with that 
observed locally. In Section 6 we present our final conclu- 
sions. Throughout this work a cosmology of Ho = 70 km 
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s 1 Mpc \ f2 TO = 0.3 and Qa = 0.7 is assumed and all mag- 
nitudes are quoted in the AB system (Oke & Gunn 1983). 



2 THE DATA 

All of the galaxies studied in this work have been selected 
from the UKIRT Ultra Deep Survey (UDS). The UDS is 
the deepest of five surveys being undertaken with the UK 
Infra-Red Telescope in Hawaii, which together comprise the 
UK I nfrared Deep Sky Survey (UKIDSS; 
2007}. The UDS covers an area of 0.8 sq 
Dec= 



Lawrence et al. 



degrees, is cen- 

tred on RA=02:17:48, Dec=-05:05:57 and is currently the 
deepest, large area, near-IR survey in existence. In addition 
to the available UKIRT near-IR imaging, the UDS is sup- 
ported by a wide array of multi-wavelength datasets, ranging 
from the X-ray to radio regimes (e.g. Cirasuolo et al. 2010). 
In this section we briefly summarize the properties of the 
UDS datasets which have been directly exploited during the 
course of this study. 

2.1 Photometric data 

The latest release of the UDS, data release eight (DR8), con- 
sists of JHK imaging to 5a— depths of J = 24.9, H = 24.3 
and K — 24.7 (2 arcsec-diameter apertures). Moreover, in 
this study we have also made use of Y— band imaging of the 
UDS (Y=23.6; 5a) obtained during the verification phase 
of the VISTA VIDEO survey (Jarvis et al. 2012, in prepa- 
ration). In addition to the available near-IR imaging, the 
second key imaging dataset in the UDS is the deep Sub- 
aru optical imaging obtained as part of the Subaru/XMM- 
Newton Deep Survey (SXDS). The SXDS imaging (Furu- 
sawa et al. 2008) covers a total area of 1.3 sq. degrees to 
5a depths of B = 27.6, V = 27.2, R = 27.0, i' = 27.0 
and z' = 26.0 (2 arcsec-diameter apertures). The final key 
dataset used in this work is the Spitzer Public Legacy Survey 
of the UKIDSS Ultra Deep Survey (SpUDS). The SpUDS 
imaging covers the full near-IR survey area and consists of 
IRAC imaging to 5a depths of 3.6/im=23.7, 4.5/im=23.3, 
5.8/im=21.6, 8.0/im=21.4 (3.8 arcsec-diameter apertures) 
and MIPS imaging to a 5a depth of 24/xm=18.9 (15 arcsec- 
diameter apertures). 

2.2 Spectroscopic data 

The optical spectra analysed here were observed on the VLT 
as part of the systematic spectroscopic follow-up of the UDS 
obtained through the ESO large programme ESO 180. A- 
0776 (UDSz; P.I. O. Almaini). Although full details of UDSz 
will be presented in Almaini et al. (2012, in preparation), 
below we provide a brief description of the key elements 
within the context of this study. 

The UDSz programme was allocated a total of 235 hours 
of observations, with 93 hours allocated for observations 
with the VIMOS spectrograph (8 pointings) and 142 hours 
allocated for FORS2 observations (20 pointings). The pri- 
mary science driver for UDSz was to obtain spectroscopic 
observations of a representative sample of if— band selected 
galaxies (K < 23), photometrically pre-selected to lie at 
redshift z p hot ^ 1- Consequently, the selection of primary 
targets was designed to obtain spectra for a random l-in-6 



sampling of K < 23 galaxies with z p h a t ^ 1 (with a control 
sample of ~ 500 galaxies with z p hot < 1). 

In order to exploit the different strengths of the two 
spectrographs, bluer galaxies (V ^ 25) were targeted with 
VIMOS and redder galaxies (V > 25 A i < 24.5) were tar- 
geted with FORS2. All of the galaxies analysed in this pa- 
per were part of the spectroscopic sample observed with 
FORS2 and each target received 5 hours of on-source 
integration with the medium-resolution GRS-300I grism 
(6000A < A < 10000A with R = 660). In total 718 
(K < 23) galaxies were targeted with FORS2, returning 
451 science-grade redshifts. 



2.3 The K-Bright galaxy sample 

The sample analysed in this study consists of all eighty- 
one UDSz galaxies with robust FORS2 spectroscopic red- 
shifts (i.e. based on multiple spectral features) in the range 
1.3 ^ z ^ 1.5 which are brighter than Ktot 21.5 (hereafter 
the K-Bright sample). This magnitude limit is bright enough 
to fully exploit the large-area coverage of the UDS and to 
ensure that each member of the sample has an unambigu- 
ous spectroscopic redshift, while still being deep enough to 
produce a final sample which is mass-complete to a limit of 
M* = 6 x 10 10 Mq, even for the most passive galaxies (see 
Section 4). The upper redshift limit of z = 1.5 was chosen 
to guarantee coverage of the 4000A break at the red end of 
the FORS2 spectra and the lower redshift limit of z = 1.3 
was adopted to ensure that the FORS2 spectra also cover 
the age-sensitive spectral breaks at 2600A and 2800A. 



2.3.1 A representative sample 

If the results derived from this study are to be applied to the 
general massive galaxy population at z ~ 1.4, it is clearly 
necessary to demonstrate that the K-Bright sample is not 
a highly biased sub-set of the general galaxy population at 
this epoch. In order to investigate this issue, in the upper 
panel of Fig. 1 we plot the distribution of i — K colours for 
the K-Bright sample, compared to that of all N=1178 galax- 
ies in the UDS photometric redshift catalogue of Cirasuolo 
et al. (2010) which satisfy the criteria: 1.3 ^ z p hot ^ 1.5 and 
Ktot ^ 21.5. Although similar, a two-sample KS test demon- 
strates that the two distributions are actually statistically 
distinguishable at the ^ 3a level. As can clearly be seen 
from Fig. 1, unusually for a spectroscopic sample, the differ- 
ence is caused by the fact that the K-Bright sample does not 
feature a large enough tail of blue (i — K < 2) objects. This 
feature of the FORS2 spectroscopic sample is by design and 
is the result of deliberately targeting the bluer members of 
the UDSz sample with the VIMOS spectrograph. 

In the lower panel of Fig. 1 we show the i — K distri- 
bution of the K-Bright sample compared to the N=970 ob- 
jects from the Cirasuolo et al. photometric redshift catalogue 
which satisfy: 1.3 z p hot ^ 1.5, K to t 21.5 andi — K ^ 2.2. 
With this additional colour criterion the two distributions 
are now statistically indistinguishable (p = 0.22). For galax- 
ies at z ~ 1.4 with Ktot 21.5, insisting on a colour of 
i — K ^ 2.2 is effectively the same as insisting that stellar 
mass is M t ^ 3x 10 10 Mq (see Section 4). Consequently, we 
conclude that the K-Bright sample is a statistically random 
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2.3.2 Sampling factor 

The UDSz spectroscopic programme was designed to pro- 
vide a random, l-in-6 (/ romp = 6), sampling of the z 1 
galaxy population with K ^ 23. However, at the Ktot ^ 21.5 
magnitude limit of the K-Bright sample the FORS2 compo- 
nent of the UDSz spectroscopy programme actually provides 
a significantly higher sampling rate of the very high-mass 
end of the galaxy mass function at 1.3 < z < 1.5. Based 
on the Cirasuolo et al. photometric redshift catalogue, and 
accounting for the differences in relative areal coverage, we 
calculate that, compared to the underlying galaxy popula- 
tion with 1.3 < z v hot < 1.5, K tot ^ 21.5 and i — K ^ 2.2, the 
sample factor of the K-Bright sample is f sa mp = 3.7±0.4. It 
is this sampling factor which is adopted when comparing the 
K-Bright sample to the local galaxy stellar mass function in 
Section 5. 



i-K 



3 GALAXY MASSES, AGES AND SIZES 




i-K 

Figure 1. The top panel shows the (normalised) distribution 
of the i — K colours of the eighty-one galaxies in the K-Bright 
sample (solid red line) over-plotted on the i — K distribution 
for all N=1178 galaxies in the UDS with Ktot < 21.5 and 
1.3 

z phot ^ 1-5. The two distributions are statistically distin- 
guishable at the ^ 3a level based on a two-sample KS test, which 
is due to the fact that, by design, the FORS2 sample does not 
feature the tail of blue objects with i — K < 2. The bottom panel 
compares the i — K distribution of the K-Bright sample to the 
N=970 galaxies in the UDS with K tot < 21.5, 1.3 < z phot < 1.5 
and i — K ^ 2.2. The two distributions shown in the bottom panel 
are statistically indistinguishable (p = 0.22). Given the tight re- 
lationship between i — K colour and stellar mass for galaxies at 
z ~ 1.4 with Ktot ^ 21.5, the bottom panel demonstrates that 
the K-Bright sample is a random sample of the general z ~ 1.4 
galaxy population at stellar masses Mi, 5s 3 X 1O 1O M . 



In addition to a precise spectroscopic redshift, for each mem- 
ber of the K-Bright sample there is a deep, flux calibrated, 
red-optical FORS2 spectrum and excellent multi- wavelength 
photometry spanning the range 0.4/mi < A < 4.5/im. Con- 
sequently, in order to take full advantage of the available 
information, we have derived stellar mass and age measure- 
ments for each member of the K-Bright sample via simul- 
taneous spectro-photometric fitting of the FORS2 spectrum 
and multi-wavelength photometry. 

3.1 Spectro-photometric fitting 

For the purposes of the spectro-photometric fitting we adopt 



both the |Bruzual fc Charlot1j2003[ ) (hereafter BC03) and the 
Chariot & Bruzual (2007^ (hereafter CB07) stellar pop- 
ulation synthesis models. The CB07 models were chosen 
because they include an updated treatment of thermally 
pulsating asymptotic giant branch (TP-AGB) stars which 
could, at least potentially, have a large influence at rest- 
frame near-IR wavelengths for stellar populations with ages 
of 0.3 ^ t ^ 2 Gyr (e.g. Maraston 2005). However, given 
that the influence of TP-AGB stars remains controversial 
(e.g. Kriek et al. 2010), it was decided to investigate the 
range of stellar masses and ages returned by both the BC03 
and CB07 stellar population synthesis models. 

Moreover, as galaxies are complex systems, potentially 
consisting of several stellar sub-systems from different gen- 
erations of star- formation, in some cases it is necessary to 
invoke models with a variety of star-formation histories in 
order to obtain good fits to the data. To account for this, 
the data for each member of the K-Bright sample was fit- 
ted using an extensive set of models encompassing a range 
of different star-formation histories, metallicities and dust 
attenuations. 



sub-sample of the galaxy population at 1.3 ^ z ^ 1.5 with 
stellar masses ^ 3 x 10 10 M Q . For information, in Table 
Al of the appendix we provide the basic observational data 
and derived properties for the eighty-one members of the 
K-Bright galaxy sample. 



3.1.1 Exponentially decaying star-formation rate models 

The data for each member of the K-Bright sample was 
initially fitted using BC03 and CB07 models with ex- 

1 Stephane Chariot, private communication 
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Figure 2. Plots illustrating the spectro-photometric fitting process for four example objects from the K-Bright sample. The left-hand 
panels show the combined fits to the spectra and photometry, while the right-hand panels show the model fits through the spectra in 
detail. In each case the optical spectrum is shown in black, the error spectrum is shown in green, the photometry is in red and the 
best-fitting model is in grey. For presentation purposes the optical spectra have been smoothed using a 5-pixel boxcar. 



ponentially decaying star-formation rates (SFR); so-called 
"t— models". The r— models are parameterized such that 
SFR oc exp(— t/r), where t is the time since star forma- 
tion began and r is the e— folding time. A wide range of 
e— folding times were considered, with the extreme values 
of r = 10 Gyr and t — 0.1 Gyr approximating constant 
star-formation and instantaneous burst models respectively. 
As listed in Table 1, in addition to e— folding time, a wide 
range of ages, metallicities and dust attenuations were also 
explored. All models employed in the spectro-photometric 



fitting were based on a Chabrier (2003) IMF and assumed 



the Calzetti et al. (2000) dust attenuation curve. 



3.1.2 Double burst models 

Although the r— models effectively describe a range of differ- 
ent star-formation histories, they do not adequately describe 
the situation whereby the majority of a galaxy's stellar mass 
is produced rapidly at high redshift (and dominates rest- 
frame near-IR wavelengths) but a low-level of more recent 
star-formation dominates the observed flux at rest-frame 
UV/optical wavelengths. In order to account for this pos- 
sibility, the data for each member of the K-Bright sample 
was fitted using so-called "double burst models" (hereafter 



DB models). The DB models consist of two sequential in- 
stantaneous bursts (simple stellar populations), with one 
component required to be "old" (age 55 0.5 Gyrs) and one 
component required to be "young" (age < 0.5 Gyrs). 



The DB models were constructed by normalising the 
flux of the younger stellar population to match that of the 
older population at a rest-frame wavelength of 5000A. In this 
way, when the components were combined, the two popula- 
tions naturally accounted for different fractions of the total 
galaxy mass. The mass fraction of the young component 
was allowed to vary between and 1, with a fine grid of 
values adopted at low-mass fractions due to the substan- 
tial impact on the rest-frame UV/optical from even small 
amounts of recent star-formation. Finally, dust reddening 
was applied to the combined young+old stellar population 
using the Calzetti (2000) dust attenuation law. In addition 
to providing improved flexibility, the DB models offer the 
advantage of providing an estimate of the likely maximal 
stellar mass and age for each member of the K-Bright sam- 
ple. 
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Figure 3. An illustration of the different stellar-mass distributions for the K-Bright galaxy sample obtained via spectro-photometric 
fitting based on different stellar population models and star-formation histories. The stellar mass distributions are ordered (left-to-right) 
by increasing median stellar mass (vertical solid lines) and correspond to r— models and double-burst models based on the CB07 stellar 
population models followed by the equivalent models based on the BC03 stellar populations models. From left-to-right the median stellar 
mass of the sample for each distribution are as follows: 6.3 X 10 10 M , 7.4 X 10 10 M , 8.5 X 10 10 M and 9.3 X 10 10 Mq. For the purposes 
of the analysis performed in this paper the stellar-mass measurements derived from the BC03 r— models have been adopted. However, as 
illustrated here, the difference in median mass between the extreme combinations of template library/star-formation history only differ 
at the 0.17 dex level, insufficient to affect any of the science results derived in the latter sections of the paper. 



Table 1. The parameter space of stellar population models used 
in the fitting of the K-Bright sample, where tu corresponds to 
the age of the Universe at a given redshift. 



r— model parameters: 

e— folding time: 
Age: 

Metallicity: 
Reddening: 



0.1 Gyr < t < 10 Gyr 
0.01 Gyr <t<t H 
0.2Z Q < Z < 2.5Z 
0.0 < A v < 3.0 mags 



DB model parameters: 

Old comp metallicity: 

Old comp age: 

Young comp metallicity: 

Young comp age: 

Young comp mass fraction: 

Reddening: 



O.2Z ,Z 
0.5 Gyr < toe < t H 
O.2Z < Z < 2.5Z 
0.01 Gyr < t yc < 0.5 Gyr 
0.0 < m yc < 1.0 
0.0 < A v < 3.0 mags 



3.1.3 Model fitting procedure 

For each member of the K-Bright sample, the stellar popula- 
tion models were initially shifted to the appropriate redshift 
and the corresponding FORS2 spectrum was binned-up to 
match the resolution of the stellar population model (typ- 
ically 10 A per pixel). The stellar population models were 
then simultaneously fitted to the spectra and the accompa- 
nying BVYJHK and IRAC 3.6/^m & 4.5/j.m photometry. 
During the fitting process the FORS2 spectral data in the 
wavelength range 6200-9600A were included in the fit, but 
the wavelength range 9600-10000A was excluded to ensure 
that bad sky-line residuals did not dominate the fit. The 
photometric data covering the same wavelength range as the 
spectra (Riz) was not included in the fitting process, simply 
because those photometric bands had already been used to 
flux normalize the spectra and photometric datasets. 

The photometric broad-band magnitudes were con- 
verted to flux densities at the effective wavelength of each 
filter and the corresponding model predictions were calcu- 



lated by integrating the models over the appropriate filter 
profiles. Based on the parameter space outlined in Table 1, 
each galaxy was fitted using a grid of up to ~ 5 x 10 mod- 
els from both the BC03 and CB07 libraries. The best-fitting 
model was determined using x 2 minimization, with the pho- 
tometry weighted using the appropriate photometric errors 
and the spectral data weighted using the FORS2 error spec- 
tra. 



In order to prevent the fitting procedure from being 
dominated by the spectral data two key steps were neces- 
sary. Firstly, it was necessary to mask areas of the spectra 
which the spectral synthesis models were not capable of re- 
producing (principally the [Oil] emission line). Secondly, it 
was necessary to ensure that the FORS2 error spectra ac- 
tually accounted for all the inherent uncertainties present 
in the data (both random and systematic). The final re- 
duced FORS2 spectrum for each object typically consists of 
5 hours of on-source integration time, comprised of a stack of 
twenty-four individual exposures, each of 12.5 minutes inte- 
gration. Consequently, at each spectral pixel it was possible 
to use these independent flux measurements to calculate the 
sigma (of) which fully accounted for the variation in mea- 
sured flux due to shot noise, varying observing conditions 
and systematic errors such as imperfect sky-line subtrac- 
tion. By weighting the final FORS2 spectra by the standard 
error (i.e. <7f/v24), it was possible to perform a combined 
X 2 fit to the spectra+photometry without either component 
dominating the fit. 



Four examples of the best-fitting model spectral energy 
distributions (SEDs) are shown in Fig. [5] where the optical 
galaxy spectra are plotted in black, the error spectra are 
plotted in green, the photometric data-points are plotted in 
red and the best-fitting models are plotted in grey. For each 
object, the left-hand panel of each row shows the full fit and 
the right-hand panel shows a close-up of the fitting region 
containing the optical spectra. For presentation purposes the 
optical spectra have been smoothed using a 5-pixel boxcar. 
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5. 1.4 Adopted stellar masses 

In Fig. 3 we show the stellar-mass distributions for the K- 
Bright sample, as derived using (left-to-right) t— models and 
DB models based on the CB07 template libraries, followed 
by the corresponding r— models and DB models based on the 
BC03 template libraries. Although, as expected, the models 
based on the CB07 templates produce lower stellar-mass es- 
timates than the models based on the BC03 templates (and 
the r— models produce lower stellar-mass estimates than the 
DB models), perhaps the most striking feature of Fig. 3 is 
the similarity between the resulting stellar-mass distribu- 
tions. This is highlighted by the fact that even between the 
two extremes, the median stellar mass of the K-Bright sam- 
ple only changes by 0.17 dex. It is clear from Fig. 3 that 
the TP-AGB phase is not having a significant impact on 
the derived stellar masses of the K-Bright sample. Further- 
more, the similarity of the derived stellar-mass distributions 
implies that the proceeding exploration of the size-mass re- 
lation at z ~ 1.4 will not be dominated by systematic stellar- 
mass uncertainties driven by galaxy SED template choice. 
Armed with this information, throughout the rest of the pa- 
per we simply adopt the stellar-mass and age measurements 
derived from the BC03 t— models, in order that they can be 
directly compared to numerous other studies in the litera- 
ture (see Section 4.5). 
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Figure 4. A comparison between the effective half-light radii 
derived from ground-based K— band imaging and HST H— band 
imaging for a sample of N=22 objects with Ktot ^ 21.5 and 
1.3 < z phot < 1.5 within the CANDELS/UDS field. The thick 
red line shows the best-fitting relationship (Equation 1) and the 
thin dashed line is a one-to-one relationship (see text for a full 
discussion) . 



3.2 Galaxy size measurements 

Measurements of the galaxy radii were derived using two- 
dimensional modelling of the latest (DR8) K— band imag- 
ing of the UDS which has a 5a detection limit of K = 24.7 
(2" diameter apertures) and 0.75" FWHM image quality. 
Over the redshift range covered by the K-Bright sample 
(1.3 ^ 2 ^ 1.5) the UDS K— band imaging is centred on 
rest-frame wavelengths of 8800A < A < 9600A, which is 
advantageous for two reasons. Firstly, at these rest-frame 
wavelengths the K— band imaging provides a size measure- 
ment significantly long-ward of the 4000A break, where the 
galaxy SED should be relatively unaffected by any recent 
episodes of star-formation. Secondly, at 2 ~ 1.4 the UDS 
K—band imaging provides size measurements well matched 
in terms of rest-frame wavelength to the benchmark galaxy 
size measurements from Shen et al. (2003), which are based 
on z— band SDSS imaging. 

The galaxy radii were calculated using a modified ver- 



sion of GALAPAGOS (Barden et al 2012) and GALFIT (Peng 



et al. 2002). GALAPAGOS is a wrapper script for GALFIT which 
fits single Sersic profiles to each object in the image, taking 
initial parameters from a catalogue generated by SExtrac- 
tor (Bertin & Arnouts 1996). A key advantage of GALAPA- 
GOS compared to GALFIT alone is that it has an algorithm 
for background estimation and, if necessary, will simulta- 
neously model any nearby companion objects which would 
otherwise unduly influence the fit. 

For this study the code was adapted to produce a 
densely sampled grid of R e , Sersic index and background 
level in order to estimate accurate errors on the fitted pa- 
rameters from contours of constant A^ 2 - It is worth empha- 
sising that, because both R e and Sersic index correlate with 
the sky-background level, including a variable background 



in the parameter estimation grid is essential for obtaining 
realistic uncertainties on the derived parameters. 



3.2.1 Comparison with higher resolution imaging 

One approach to testing the accuracy of our galaxy size mea- 
surements would be to analyse simulated images containing 
artificial galaxies spanning a range of half-light radii and 
Sersic index. We have not performed such a simulation in 
this work, simply because building and then reclaiming the 
properties of artificial galaxies with the same axi-symmetric 
galaxy models is likely to over-estimate the accuracy with 
which the galaxy sizes can be recovered. 

However, within the UDS field we are in the fortu- 
nate position that it is possible to directly compare galaxy 
sizes derived from ground-based and HST imaging in the 
near-IR. Consequently, in order to test whether or not our 
determinations of the galaxy radii have been affected by 
the ground-based spatial resolution of the UDS K— band 
imaging, we have made use of the publicly available Cos- 
mic Assembly Near-infrared Deep Extragalactic Legacy Sur- 
vey (CANDELS, Grogin et al. 2011; Koekemoer et al. 2011) 
imaging in the UDS. The CANDELS dataset in the UDS fea- 
tures J- and #-band WFC3/IR imaging (FWHM ~ 0.2") 
over an area of ~ 200 sq. arcmins to a 5a— depth of 26.8 
(0.6"diameter apertures). 

Based on the photometric redshift catalogue of Bruce 
et al. (2012), there are twenty-two galaxies within the CAN- 
DELS/UDS area with Ktot ^ 21.5 and photometric redshifts 
in the range 1.3 ^ 2 ^ 1.5. For this sample the half-light 
radii and Sersic index measurements derived by Bruce et 
al. (2012), based on a GALFIT analysis of the H— band HST 
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imaging, were compared to those derived using the ground- 
based K— band imaging. 

As can been seen from Fig. 4, the results of this compari- 
son demonstrate that the two sets of half-light radii measure- 
ments are well correlated (r s = 0.95) and follow an essen- 
tially one-to-one relation. The best-fitting relation between 
the two measurements of half-light radius was found to be: 

i?wPC3 = (0.96 ± O.OS)^'^ - 08 (1) 

where i?wFC3 and Ruus are the circularised half-light radii 
(in kpc) derived from the WFC3/IR and ground-based data 
respectively. Moreover, the corresponding Sersic index mea- 
surements were also found to be well correlated (r s = 0.82), 
with only 1/22 objects returning a significantly different 
best-fitting Sersic index (hwfc = 2.0 versus tiuds = 3.5). 
Consequently, for the purposes of the proceeding analysis 
the measurements of half-light radii and Sersic indices de- 
rived from the ground-based K— band imaging are adopted 
without correction. 



4 THE STELLAR MASS-SIZE RELATION 

In this section we present the basic results of our study of the 
galaxy size-mass relation at z ~ 1.4, exploring the location 
of the K-Bright sample on the size-mass plane as a func- 
tion of specific star-formation rate, morphology and stellar 
population age. In what follows we will consistently com- 
pare the results for the K-Bright sample with the size-mass 
relations derived by Shen et al. (2003) for early (n ^ 2.5) 
and late-type (n < 2.5) galaxies in the SDSS. In order to 
perform this comparison we have converted the semi-major 
axis half-light radii (R e ) derived by GALFIT to their equiva- 
lent circularized values (R c ) using the standard conversion: 
R c = \J (b/a)R e , where a/b is the best-fitting axial ratio. 

In the left-hand panel of Fig. 5 we show the full K- 
Bright sample on the size-mass plane. Also shown are the 
size-mass relations for late-type (thin blue line) and early- 
type (thick red line) galaxies in the SDSS as derived by Shen 
et al. (2003), with the dotted lines indicating the la scatter 
associated with both relations. The shaded region at M* > 
6 x 10 10 Mq highlights the regime where the K-Bright sam- 
ple is fully mass complete. This mass completeness threshold 
has been calculated in a deliberately conservative manner 
and is based on a model galaxy at z = 1.4 with Ktot = 21.5, 
zero current star formation and an age of 4 Gyr. Conse- 
quently, at Mi, > 6 x 10 10 Mq the K-Bright sample should 
be mass complete, even for entirely passive galaxies. In the 
stellar mass interval 3 x 10 10 M Q < M* < 6 x 10 10 M the 
K-bright sample may be partially incomplete for the most 
passive galaxies, but should still be complete for actively 
star-forming galaxies, unless they are heavily reddened. The 
apparent cut-off at stellar masses of ~ 3 x 10 10 M is 
thus a function of the adopted K to t = 21.5 magnitude limit 
and the red colour selection of the FORS2 spectroscopic 
sample (i.e. i — K ^ 2). 

4.1 Specific star- formation rate 

In the right-hand panel of Fig. 5 the K-Bright sample has 
been split on the basis of specific star- formation rate (sSFR), 
where objects with sSFR < 0.1 Gyr -1 are plotted as red 



filled circles and objects with sSFR > 0.1 Gyr -1 are plotted 
as open blue circles. Although it is obviously possible to 
derive a value of sSFR directly from the best-fitting SED 
templates, in order to avoid too much coupling between the 
adopted stellar mass and star-formation rate measurements, 
the decision was taken to exploit the different empirical star- 
formation indicators available for the K-Bright sample. 

The first star-formation indicator we have employed is 
the UV luminosity calibration of Madau et al. (1998) and 
was based on the observed B— band magnitudes, dust cor- 
rected using the best-fitting value of A v . The second indi- 
cator is based on the [Oil] emission line fluxes (or upper 
limits), and is based on the Kennicutt (1998) calibration 
assuming that the Ha emission line undergoes an average 
extinction of one magnitude. The final star-formation indi- 
cator is based on the observed 24/im fluxes (or upper limits), 
which are initially converted into rest-frame 8/im luminosi- 
ties using a 10 11 Lq starburst template from Lagache et al. 
(2003). Based on the 8/um luminosity the total infra-red lu- 
minosity is then estimated using the calibration of Bavouzet 
et al. (2008) and the star-formation rate is estimated using 
the calibration of Kennicutt (1998). 

The adopted star-formation rate estimate for each mem- 
ber of the K-Bright sample is taken as the largest value de- 
rived from [Oil] emission-line flux, dust-corrected UV lumi- 
nosity and MIPS 24/im flux. In this fashion we have delib- 
erately adopted the maximum likely value of sSFR for each 
object, which should ensure that objects with extremely low 
values of sSFR can be regarded as genuinely passive. The 
choice of 0.1 Gyr -1 as the threshold between actively star- 
forming and passive objects is motivated by several different 
factors. Fundamentally, this value is the median sSFR value 
of the K-Bright sample and is also a suitable dividing line 
separating galaxies which are on or off the so-called "main 
sequence of star-formation" at this epoch (e.g. Elbaz et al. 
2011). Moreover, at z ~ 1.4 a sSFR value of 0.1 Gyr" 1 im- 
plies that a galaxy will require a Hubble time in order to 
double its stellar mass. 

4-1.1 Passive galaxies 

It can clearly be seen from the right-hand panel of Fig. 5 that 
the passive members of the K-Bright sample follow a distinct 
size-mass relation from the actively star-forming members. 
This point is clarified further in Fig. 6 which shows just 
the passive members of the K-Bright sample on the size- 
mass plane. It is apparent from the left-hand panel of Fig. 
6 that the passive members of the K-Bright sample follow a 
well defined size-mass relation, with a slope consistent with 
the size-mass relation displayed by local early-type galaxies, 
but with systematically lower half-light radii than their local 
counterparts. 

In this context it is instructive to calculate the size 
growth factor (f g ), defined here as the increase in half-light 
radius required to place the K-Bright galaxies onto the ap- 
propriate local galaxy size-mass relation from Shen et al. 
(2003). To reconcile the passive members of the K-bright 
sample with the local early-type size-mass relation requires 
a median growth factor of f g = 2.43 ± 0.20, even allowing 
for no increase in stellar mass. This off-set is illustrated by 
the right-hand panel of Fig. 6, where the local size-mass rela- 
tion for early-type galaxies derived by Shen et al. (2003) has 
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Figure 5. The galaxy size - stellar mass relation for the K-Bright sample. The left-hand panel shows the location of the K-Bright 
sample on the half-light radius - stellar mass plane. The thin (blue) and thick (red) solid lines show the relations derived for z < 0.1 
SDSS galaxies by Shen et al. (2003) for late-type (n < 2.5) and early-type (n ^ 2.5) galaxies respectively (dotted lines indicate the 
associated la scatter). The shaded region at log(M*/MQ) ^ 10.8 highlights where the K-Bright sample should be 100% mass-complete, 
even for the most passive galaxies (see text for discussion). The right-hand panel shows the same information, except that the K-Bright 
galaxies have been split into passive (red solid circles) and star-forming (blue open circles) sub-samples depending on whether their 
specific star-formation rate (sSFR) lies above or below 0.1 Gyr - 1 . In both panels, and all subsequent figures, the plotted half-light radii 
measurements have been circularized (i.e. R c = \J (b/a)R e ) in order to facilitate direct comparison with the results of Shen et al. (2003). 
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Figure 6. The left-hand panel shows the galaxy size - stellar mass relation for those members of the K-Bright sample classified as 
passive according to our adopted sSFR ^0.1 Gyr -1 threshold. The solid and dotted lines show the size-mass relation for early-type 
galaxies in the SDSS derived by Shen et al. (2003) and the corresponding la scatter respectively. It can be seen that the passive K-Bright 
galaxies follow a size-mass relation which has an identical slope (and scatter) to that of local early-type galaxies, but is simply shifted 
in normalisation. This is illustrated by the right-hand panel which shows the effect of reducing the vertical normalization of the Shen et 
al. early-type relation by a growth factor of f g = 2.43 ± 0.20 (see text for discussion). 



been lowered in normalization by a factor of f g — 2.43. Re- 
markably, it can be seen that this simple renormalization of 
the local early-type size-mass relation provides an excellent 
description of the size-mass relation followed by the passive 
members of the K-bright sample at z ~ 1.4, both in terms of 
slope and scatter. We will return to this point in Section 5. 



4-1-2 Star-forming galaxies 

In the left-hand panel of Fig. 7 we plot the star-forming 
(sSFR > 0.1 Gyr -1 ) half of the K-Bright sample on the size- 
mass plane and compare with the local size-mass relation for 
late- type galaxies in the SDSS derived by Shen et al. (2003). 
In contrast to the passive members of the K-Bright sample, 
it can immediately be seen that ~ 50% of the star-forming 
galaxies are already consistent with the size-mass relation 
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Figure 7. The left-hand panel shows the galaxy size - stellar mass relation for the those members of the K-Bright sample classified as 
star-forming according to our adopted sSFR > 0.1 Gyr - 1 threshold. The solid and dotted lines show the size-mass relation for late-type 
galaxies in the SDSS derived by Shen et al. (2003) and the corresponding lcr scatter respectively. It can be seen that the star-forming 
K-Bright galaxies are consistent with following a size-mass relation similar to that of local late-types, but shifted in normalisation. This 
is illustrated by the right-hand panel which shows the effect of reducing the vertical normalization of the Shen et al. late-type relation 
by a growth factor of f g = 1.61 ± 0.17 (see text for discussion). 



of local late-type galaxies. As a result, the median value of 
the growth factor for the star-forming members of the K- 
Bright sample is only f g — 1.61 ± 0.17^] as illustrated in 
the right-hand panel of Fig. 7. 

The results of this sub-section are broadly consistent 
with many previous studies in the literature (see Section 1) 
in that, at a given epoch, the off-set from the local size-mass 
relations is significantly larger for passive galaxies than for 
their star-forming contemporaries. A detailed comparison 
with relevant results from recent literature will be presented 
in Section 4.5. 



4.2 Morphology 

When considering the evolution of the galaxy size-mass re- 
lation it is conventional to consider the influence of galaxy 
morphology. Indeed, much of the discussion in the literature 
has concentrated on the size evolution of morphologically 
selected early-type galaxies. Consequently, in this section 
we investigate the size-mass relation for the K-Bright sam- 
ple, splitting the sample into late-type and early-type sub- 
samples using the same Sersic index threshold of n = 2.5 
adopted by Shen et al. (2003). 

4-2.1 Early-type galaxies 

The middle panel of Fig. 8 shows the size-mass relation for 
the early-type (n 2.5) members of the K-Bright sample 
compared to the size-mass relation for local early-type galax- 
ies from Shen et al. (2003). It can be seen that the early- 
type galaxies follow a size-mass relation which has a slope 

2 The outlier at (log(M/M©) = 10.5, log(R c /kpc) = -0.5) has 
been excluded from this calculation. 



which is consistent with the local size-mass relation for early- 
types, but is off-set to smaller radii at a given stellar mass. 
In terms of the growth factor (f g ), the early-type galaxies 
are off-set from the local size-mass relation by a factor of 
f g = 2.37 ±0.29. 

If we confine our attention to the growth factor required 
to reconcile them with the local size-mass relation for early- 
type galaxies, it is clear that the passive and early-type sub- 
samples of the K-Bright galaxies are very similar, with f g — 
2.43 ± 0.20 and f g = 2.37 ± 0.29 respectively. However, it 
is noticeable that the scatter associated with the size-mass 
relation displayed by the early-type sub-sample is somewhat 
larger than that associated with the size-mass relation of the 
passive galaxy sub-sample. Comparing Fig. 6 with Fig. 8 it 
is clear that the reason for the increased scatter is that, 
unlike the passive galaxy sub-sample, ~ 30% of the early- 
type sub-sample have half-light radii which are consistent 
with the local early-type size-mass relation. 

Interestingly, although splitting the full K-Bright sam- 
ple at n 2.5 or sSFR ^ 0.1 Gyr" 1 delivers similar sized 
sub-samples (N=37 and N=41 galaxies respectively), the 
correspondence between an early-type morphology and pas- 
sivity is not one-to-one. Indeed, the forty-one galaxies in 
the K-Bright sample classified as passive are split relatively 
evenly between early-type (N=23) and late-type morpholo- 
gies (N=18). Similarly, of the thirty-seven objects in the 
early-type sub-sample, only twenty-three are also classified 
as passive. Consequently, it is clear that a relatively large 
fraction (~ 20%) of the K-Bright galaxies are classified as 
star-forming early-types according to our adopted sSFR and 
Sersic index criteria. 

4-2.2 Late-type galaxies 

In the bottom panel of Fig. 8 we plot the size-mass relation 
for the late-type (n < 2.5) members of the K-Bright sample. 
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Figure 8. The top panel shows the size-mass relation for the K- 
Bright sample with early-type galaxies (n ^ 2.5) plotted as solid 
red circles and late-type galaxies (n < 2.5) plotted as open blue 
circles. The middle panel shows the early- type K-Bright galaxies 
compared to the local early-type size-mass relation from Shen et 
al. (2003). The bottom panel shows the equivalent comparison for 
the late-type K-Bright galaxies. 



As with the early-type sub-sample, the late-type members 
of the K-Bright sample appear to follow a size-mass relation 
which has a slope which is consistent with the equivalent 
local size-mass relation, but off-set to smaller half-light radii 
at a given stellar mass. If we calculate the off-set from the 
local size-mass relation in terms of the growth factor we 
derive a median value of f g = 2.15 ± 0.15. 

It is noteworthy that this growth factor is larger than 
the growth factor previously calculated for the star-forming 
members of the K-Bright sample (f g = 1.61 ± 0.17). Again, 
the reason for this apparent discrepancy is that the corre- 
spondence between late-type morphology and active star- 
formation is not one-to-one. In fact, of the forty- four mem- 
bers of the late- type sub-sample, twenty-six are actively star- 
forming galaxies and eighteen are galaxies which are classi- 
fied as passive according to our adopted sSFR threshold. 
The overall growth factor of f g = 2.15 ± 0.15 is therefore ef- 
fectively an average of the f g = 1.79 ±0.17 displayed by the 
star-forming late-type galaxies and the significantly larger 
growth factor of f g = 2.48 ± 0.23 displayed by the passive 
late-type galaxies. 

Therefore, not only is ~ 20% of the K-Bright sample 
comprised of star-forming bulges which are larger than their 
passive counterparts, a further ~ 20% is comprised of pas- 
sive disk-like galaxies which are smaller than their actively 
star-forming counterparts. In combination with the tight 
size-mass correlation displayed by the passive sub-sample 
in Fig. 5, this information leads naturally to the conclusion 
that the location of massive z ~ 1.4 galaxies on the size- 
mass plane is a stronger function of specific star-formation 
rate than galaxy morphology. 

Finally, it is interesting to note that ^ 40% of the pas- 
sive members of the K-Bright sample are classified as having 
a disk-like morphology (i.e. n < 2.5). This result provides 
potentially important information about the quenching of 
star-formation in massive galaxies and will be discussed in 
more detail in Section 5. 



4.3 Stellar population age 

In this section we examine the evidence for a relationship 
between galaxy size and stellar population age within the 
early- type/passive members of the K-bright sample, given 
that the confirmation of any such relationship would offer 
important constraints on competing evolutionary scenarios. 
Recent results on this issue in the literature have been con- 
troversial, with different studies finding apparently contra- 
dictory claims (e.g. Saracco et al. 2009; Trujillo et al. 2011). 
In the following discussion we have adopted the best-fitting 
age since formation (tj or ) as our age indicator, although 
quantitatively similar results are derived if the alternative 
age indicator tf or /r is employed. 

4-3.1 Passive galaxies 

The median age of the passive galaxy sub-sample is 3.3 ± 
0.2 Gyr, which is consistent with a formation redshift of 
zf or ~ 4.5. If we divide the sample into "old" and "young" 
based on this median age, then the passive sub-sample also 
cleanly separates in terms of half-light radius and stellar 
mass. For example, the median stellar mass of the "old" 
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passive galaxies is M* = (1.64 ± 0.29) x lO n M , whereas 
the median stellar mass of the "young" passive galaxies is 

= (0.79 ± 0.21) x 10 11 M Q . The corresponding median 
values for the half-light radii are 2.52 ± 0.35 kpc for the 
"old" passive sub-sample and 1.55±0.22 kpc for the "young" 
passive sub-sample. 

Several aspects of these results are noteworthy. Firstly, 
the fact that the oldest passive members of the K-Bright 
sample are systematically more massive than the youngest 
passive members indicates that the process of "downsizing" 
is already in place for massive galaxies at z ~ 1.4. Moreover, 
the difference in median half-light radius for the "old" and 
"young" passive galaxies is entirely consistent with the cor- 
relation between half-light radius and stellar mass seen for 
early-type galaxies in the SDSS by Shen et al. (2003), i.e. 
Thi oc M 0,56 (as expected given Fig. 6). 

Therefore, although we do find a factor of ~ 1.6 differ- 
ence in the median half-light radii of the "old" and "young" 
passive galaxy sub-samples, this is simply a consequence of 
the correlations between size-mass and mass-age. Indeed, 
because the passive galaxy sample follows a size-mass rela- 
tion with an identical slope to that seen at low-redshift (see 
Fig. 6), the median growth factors of the "old" and "young" 
members of the passive galaxy sub-sample are perfectly con- 
sistent (f g = 2.39 ± 0.32 and f g = 2.41 ± 0.26 respectively). 
Based on this evidence we conclude that there is no correla- 
tion between stellar population age and off-set from the local 
size-mass relation for passive galaxies within the K-Bright 
sample (see Section 5 for further discussion). 

4-3.2 Early-type galaxies 

A more direct comparison with existing studies in the liter- 
ature is to examine if there are any correlations with stellar 
population age within the early-type (n ^ 2.5) members of 
the K-Bright sample. 

To explore this issue, we adopt the same simple strat- 
egy and split the early-types into "young" and "old" sub- 
samples using the median age since formation (which is 
again 3.3 Gyr). If we then calculate the median growth fac- 
tor needed to reconcile the K-Bright early-types with the 
local early-type size-mass relation we derive values of f g = 
2.50 ± 0.32 for the "old" sub-sample and f g = 2.28 ± 0.46 
for the "young" sub-sample. As for the passive members of 
the K-Bright sample, based on these results we conclude 
that there is no indication that the off-set from the local 
early-type size-mass relation is related to stellar population 
age at z ~ 1.4 (consistent with the conclusions of Cimatti, 
Nipoti & Cassata 2012). In fact, as with the passive sub- 
sample, the median half-light radius of the "old" early-type 
galaxies is actually larger than that of the "young" early- 
types (2.2 ± 0.3 kpc versus 1.6 ± 0.3 kpc), as a result of the 
correlation between stellar population age and stellar mass. 



4.4 Dynamical mass measurements 

For a thirteen object sub-sample of the K-Bright galaxies 
it was possible to extract a reliable measurement of their 
stellar-velocity dispersion. Consequently, for this small sub- 
sample it is possible to derive dynamical-mass estimates and 
thereby provide an independent measurement of the typical 
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Figure 9. Illustration of the model fitting employed to derive 
stellar-velocity disperions for thirteen galaxies with sufficiently 
high signal-to-noise (SNR ~ 8 per pixel) in the Ca H & K spectral 
region. In each panel the thick green line is the galaxy spectrum, 
the black line is the best-fitting convolved stellar spectrum and 
the thin grey line is the error spectrum. The catalogue ID number 
and best-fitting stellar-velocity dispersion are shown at the top 
of each panel. 

size evolution, free from the systematic uncertainties associ- 
ated with stellar-mass measurements based on SED fitting. 
The sub-sample of galaxies with stellar-velocity dispersion 
measurements is confined to those objects which were bright 
enough to provide the necessary signal-to-noise (i.e. SNR 
~ 8 per pixel) and to the redshift interval 1.30 ^ z ^ 1.41, 
such that the strong Ca H & K absorption features lie within 
a clean region of spectrum, clear of strong sky-line residuals. 

The velocity dispersions were measured via \ 2 mini- 
mization in pixel space, using a template library of 215 high 
signal-to-noise stellar spectra taken from STELIB Stellar Li- 
brary (Le Borgne et al. 2003), as illustrated in Fig. 9. The 
corresponding dynamical mass estimates were derived using 
the following expression: 

M dyn = ^ (2) 

where R c is the circularised half-light radius and a is the 
line-of-sight velocity dispersion within the half-light radius 
(aperture corrected according to Jorgensen et al. 1996) . 

4-4-1 Size evolution 

Based on the relationship between half-light radius and dy- 
namical mass derived for early-type SDSS galaxies by van 
der Wei et al. (2008), it is possible to predict the z — 
half-light radius of a galaxy with a given dynamical mass. 
Therefore, for each K-Bright galaxy with a measured ve- 
locity dispersion it is possible to calculate the ratio of the 
observed half-light radius at z ~ 1.4 to the expected half- 
light radius at z — for a galaxy of the same dynamical 
mass. 

In Fig. 10 we show the results of this calculation for the 
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Figure 10. The evolution of circularised effective radius as a 
function of redshift for early-type galaxies with dynamical mass 
estimates of Mdyn — 10 11 Mq, normalized to the dynamical mass 
- size relation determined from the SDSS by van der Wei et al. 
(2008); see text for details. The orange data points are taken from 
the sample of 50 early-type galaxies published by van der Wei et 
al. (2008), which contains both field galaxies at z ~ 1 and cluster 
galaxies at z = 0.831. The green data points are a sample of 17 
early-type field galaxies in the redshift interval 1.054 ^ z ^ 1.598 
with velocity dispersion measurements derived from Keck spec- 
troscopy published by Newman et al. (2010). The black data 
points are the 13 galaxies from the K-Bright sample for which 
it was possible to obtain velocity dispersion measurements and 
therefore dynamical mass estimates. In each case the filled trian- 
gles with error bars indicate the sample medians. The solid curve 
shows the evolutionary trend derived by van der Wei et al. (2008) 
which has the functional form R c (z) <x (l + z) — 0.98±0.n_ j n orc j er 
to agree with the velocity dispersion normalization adopted both 
here and in van der Wei et al. (2008), the velocity dispersion mea- 
surements of Newman et al. (2010) have been corrected to match 
the velocity dispersion within the half-light radius according to 
the prescription outlined in Jorgensen et al. (1996). 



thirteen K-Bright galaxies with velocity-dispersion measure- 
ments. In Fig. 10 we also show the van der Wei et al. (2008) 
sample of 50 early- type galaxies, consisting of both z ~ 1 
field galaxies and cluster galaxies at z = 0.831, and the re- 
sults from Newman et al. (2010) for 17 early-type galaxies 
in the redshift interval 1.054 ^ z ^ 1.598. The solid curve 
shown in Fig. 10 is the best-fitting relation for the observed 
half-light radius evolution from van der Wei et al. (2008) 
which has the form: R c (z) oc (1 + Z )~°- 9S±0 

It can be seen from Fig. 10 that the results derived 
here are in good agreement with those derived by Van der 
Wei et al. (2008) and Newman et al. (2010). Moreover, us- 
ing the median value of R c / ' R c (z=a) for the K-bright sub- 
sample with dynamical-mass measurements we can derive 
an independent estimate of the growth factor needed to rec- 
oncile the sizes of the K-Bright galaxies with their local 
counterparts. It is encouraging that the resulting value of 
f g = 2.33 ± 0.32 is fully consistent with our measurements 
derived on the basis of the stellar mass - size relation. 



4.5 Comparison with the literature 

Given the extensive body of work dealing with the size evo- 
lution of massive galaxies already in the literature, it is 
clearly impossible to comprehensively compare our results 
with every relevant paper on the subject. Consequently, in 
this sub-section we have chosen to compare our results with 
four representative papers which cover the same redshift and 
stellar-mass range as the K-Bright sample. 



4.5.1 Trujillo et al. (2007) 

Trujillo et al. (2007) analysed the sizes of massive (M* ^ 
10 n M Q ) galaxies in the redshift interval 0.2 < z < 2.0, 
based on a sample of 810 if— band selected galaxies in 
the Palomar/DEEP-2 survey. Using stellar-mass estimates 
based on the BC03 stellar population models (Chabrier 
IMF) and circularized half-light radii derived from the avail- 
able F814W ACS imaging, Trujillo et al. observed rapid size 
evolution in both the massive late-type (n < 2.5) and early- 
type (n 2.5) galaxy populations. Given that the Trujillo 
et al. sample is if— band selected, adopts the same SED fit- 
ting models and covers the same stellar-mass range as the 
K-Bright sample, it is clearly of interest to compare their 
size evolution results with those derived here. 

A precise comparison is slightly problematic because the 
redshift range spanned by the K-Bright sample (1.3 ^ z ^ 
1.5) straddles two of the redshift bins adopted by Trujillo 
et al. (1.1 < z < 1.4 and 1.4 < Z < 1.7). However, if we 
simply adopt the average of the results for the two Trujillo 
et al. redshift bins, their results suggest that at z ~ 1.4 
the growth factor for late-type galaxies is f g ~ 1.8 and the 
growth factor for early-type galaxies is f g ~ 3.3. Exactly 
the same growth factors are obtained if we instead adopt 
the fitting formulas derived by Buitrago et al. (2008), who 
incorporated the Trujillo et al. (2007) results in their study 
of the size evolution of massive galaxies over the redshift 
interval 0.0 < z < 3.0. These figures are to be compared 
with the results presented in Section 4.2, which were f g = 
2.15±0.15 for late-types and f g — 2.37T0.29 for early-types. 

It is clear from this comparison that while our results for 
massive late- type galaxies are in reasonable agreement, our 
determination of the size evolution for massive early-type 
galaxies at z ~ 1.4 is significantly smaller than derived by 
Trujillo et al. (2007). However, in reality, the locus of the K- 
Bright early-type galaxies on the size-mass plane is entirely 
consistent with the Trujillo et al. sample (e.g. their Fig. 7). 
The fundamental reason for the difference in derived growth 
factors is simply that the K-bright sample does not display 
the significant (~ 20%) tail of objects with R c ^ 1 kpc which 
is present in the Trujillo et al. (2007) sample. It is possible 
that the lack of R c 1 kpc objects in the K-Bright sample 
is associated with the limitations of ground-based imaging, 
although our comparison to the available WFC3/IR imaging 
(see Section 3.2.1) suggests that the K— band R c determina- 
tions are reliable down to ~ 1 kpc and there is no indication 
of a "plateau" in the R c values shown in Fig. 6. 



4.5.2 Saracco et al. (2009) 

Saracco, Longhetti & Andreon (2009) investigated the mas- 
sive galaxy size-mass relation using H— band HST NIC- 
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MOS imaging of a heterogeneous sample of 32, morpholog- 
ically classified, early-type galaxies in the redshift interval 
1 < 2 < 2. Based on the available 0.4pm-2.2^im photometry, 
Saracco et al. (2009) performed SED fitting using the BC03 
template library (Chabrier IMF) and exponentially decaying 
SFHs. Consequently, it is straightforward to directly com- 
pare the results of Saracco et al. (2009) with those derived 
here. 

Based on their SED fitting, Saracco et al. found that 
their sample of early-type galaxies (ETGs) displayed a bi- 
model distribution of ages, with peaks at ~ 1 Gyr and 
~ 3.5 Gyr respectively. Splitting their sample into old ETGs 
(oETGs) and young ETGs (yETGs) using an age threshold 
of 2 Gyr, Saracco et al. found that while z ~ 1.5 oETGs 
have half-light radii which are a factor of ~ 2.5 smaller than 
their local counterparts, yETGs at z ~ 1.5 are fully consis- 
tent with the local size-mass relation. As already discussed 
in Section 4.3, amongst the early-type (n ^ 2.5) members 
of the K-Bright sample we find no evidence that the off-set 
from the local size-mass relation is a function of stellar pop- 
ulation age, instead finding that sSFR is the fundamental 
physical parameter which determines the location of mas- 
sive z ~ 1.4 galaxies on the size-mass plane. 

In light of this apparent contradiction, it is clearly of in- 
terest to ask whether we can reproduce the results of Saracco 
et al. if we employ the same criteria to separate the K-Bright 
sample. Specifically, if we identify old early-types on the ba- 
sis of n > 2.5 and a best-fitting stellar population age ^ 2 
Gyr, we find that the K-Bright sample contains twenty-nine 
such galaxies. Notably, within this sub-sample of oETGs 
22/29 (76%) would also be classified as passive according 
to our adopted sSFR threshold of < 0.1 Gyr^ 1 . Moreover, if 
we identify young early-type galaxies in the K-Bright sample 
as having n < 2.5 and age < 2 Gyr, we find that 7/8 (88%) 
of the galaxies satisfying these criteria would also be clas- 
sified as actively star-forming according to our sSFR > 0.1 
Gyr -1 criterion. 

Based on the results derived here, we would there- 
fore expect that separating the Saracco et al. (2009) sam- 
ple into oETGs and yETGs should produce a size-mass 
plot which looks similar to the right-hand panel of Fig. 
5. Specifically, we would expect that the oETGs should 
closely mimic our passive galaxy sub-sample, following a 
size-mass relation with a normalization a factor of ~ 2.5 
lower than the local early-type size-mass relation. Further- 
more, we would expect the yETGs to mimic our star-forming 
sub-sample, lying on the local early-type size-mass relation 
at log(M/M ) < 10.8, but lying a factor of ~ 2 below it 
at log(M/M0) ^ 10.8. If we examine the results derived by 
Saracco et al. (e.g. their Fig. 9) it can be seen that their sub- 
samples of yETGs and oETGs display exactly the expected 
behaviour. 

Consequently, we conclude that the Saracco et al. (2009) 
results are entirely consistent with sSFR being the funda- 
mental physical parameter governing the location of massive 
galaxies on the size-mass plane and that, amongst morpho- 
logically selected early-type galaxies at 1 < z < 2, a stellar 
population age > 2 Gyr is a reasonable, but not one-to-one, 
proxy for passivity. 



4.5.3 Williams et al. (2010) 

Williams et al. (2010) studied the size evolution of massive 
galaxies at z p hot ^ 2, using a sample of ~ 30, 000 (K < 22.4) 
galaxies selected from the first data release of the UDS. The 
stellar-mass estimates derived by Williams et al. were based 
on the BC03 template library and assume a Chabrier IMF, 
again making the comparison of their results with those de- 
rived here straightforward. 

In comparison to the Williams et al. sample, the K- 
Bright sample has the advantage of spectroscopic redshifts, 
significantly deeper photometry in the range 1.0pm— 4. 5pim 
and full spectro-photometric fitting of the combined pho- 
tometry and FORS2 spectra. However, in contrast, the 
Williams et al. sample is ~ 350 times larger than the K- 
Bright sample and covers a wider redshift range. 

As part of their investigation, Williams et al. split their 
sample into passive and star-forming galaxies using a sSFR 
threshold of sSFR= 0.3/in, where tu is the age of the Uni- 
verse at a given redshift. Notably, over the redshift range 
spanned by the K-Bright sample (1.3 ^ z < 1.5) this thresh- 
old in sSFR is virtually identical to the alternative threshold 
of sSFR= 0.1 Gyr -1 adopted here. Over the redshift interval 
0.5 < z p hot < 2.0, Williams et al. find that the evolution of 
the median half-light radii of massive (log(M/M ) > 10.8) 
passive galaxies evolves as log(i? c ) = 0.78(1 + z)~ 1A1 . 

Consequently, at the median redshift of the K-Bright 
sample, the results of Williams et al. would predict that 
the median circularized half-light radius of passive K-Bright 
galaxies with log(M/M ) > 10.8 should be 1.9 ± 0.1 kpc. 
This is in excellent agreement with the equivalent value of 
R c — 1.93 ± 0.24 kpc we derive for the passive members of 
the K-Bright sample. 

4.5.4 van Dokkum et al. (2010) 

In an effort to avoid the effects of progenitor bias, van 
Dokkum et al. (2010) used the newfirm medium band sur- 
vey (NMBS) to select galaxies with a constant number den- 
sity of n = 2 x 10~ 4 Mpc~ 3 over the photometric redshift 
interval 0.2 < z < 2.2. In each of four redshift bins, stacked 
images of the selected galaxies were used to increase the 
signal-to-noise ratio and allow the evolution of the typical 
galaxy structural parameters to be measured. Over the red- 
shift interval < z < 2.2 van Dokkum et al. (hereafter 
vDIO) found that massive galaxies selected to have a con- 
stant number density of n = 2 x 10~ 4 Mpc~ 3 display a factor 
of ~ 4 growth in their effective radius, but that the corre- 
sponding growth in their stellar mass is only a factor of 
~ 2. As noted by vDIO, these results are in agreement with 
a galaxy growth model based on minor mergers (see Sec- 
tion 5), a conclusion further supported by their finding that 
the additional stellar-mass growth occurs predominantly at 
large galactic radii (i.e. r > 5 kpc). 

Given that vDIO do not split their sample into early 
and late- type morphologies, it is interesting to compare their 
results with our results for the full K-Bright sample, in- 
dependent of morphological or sSFR thresholds. Based on 
their fitting functions for the evolution of stellar mass, vDIO 
would predict that galaxies with a constant number density 
of n — 2 x 10~ 4 Mpc~ 3 at z ~ 1.4 should have a median stel- 
lar mass of log(M/M ) = 11.24. The stellar masses for the 
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vDIO sample were calculated assuming a Kroupa IMF and 
the Maraston et al (2005) stellar population models. If we 
adopt the conversions calculated by Cimatti et al. (2008), we 
see that a stellar mass of log(M/ M ) = 11.24 is equivalent 
to log(M/M ) = 11.35 assuming a Chabrier IMF and the 
BC03 stellar population models. If we then adopt the vDIO 
policy and define our mass bin as 11.35 ± 0.15, we identify 
a sample of seventeen galaxies within the K-Bright sam- 
ple. The corresponding number density of this sub-sample 
is n — (1.5 ± 0.4) x 10~ 4 Mpc~ 3 , in good agreement with 
the vDIO expectation. However, our results indicate that 
the median effective half-light radius (semi-major axis) of 
this sub-sample is 3.3 ± 0.3 kpc, approximately twenty-five 
percent smaller than the vDIO result of 4.3 ± 0.3 kpc. 



5 DISCUSSION 

In this final section we provide an extended discussion of 
the two most important results of this study. Firstly, we 
explore whether any of the commonly proposed evolutionary 
mechanisms can plausibly explain the observed size-mass 
evolution of the passive members of the K-Bright sample, 
given the available observational and theoretical constraints. 
Secondly, we discuss the implications for the quenching of 
star-formation in massive galaxies at high-redshift given our 
result that a substantial fraction of the passive members of 
the K-Bright sample at z ~ 1.4 have disk-like morphologies. 

5.1 Mechanisms for galaxy size evolution 

Perhaps the cleanest result to emerge from this study is that 
the passive (sSFR < O.lGyr" 1 ) members of the K-Bright 
sample follow a size-mass relation which is identical in slope 
(and scatter) to that followed by early-type galaxies in the 
SDSS, but off-set to smaller half-light radii by a factor of 
~ 2.4. As a consequence, in this sub-section we investigate 
whether any of the commonly-discussed evolutionary mech- 
anisms are capable of reconciling the size-mass relations for 
passive galaxies at z ~ 1.4 and z ~ 0, without violating 
other observational constraints at low redshift. 

In the following discussion we assume that the ad- 
ditional stellar mass added by on-going star-formation 
amongst the passive galaxy sub-sample is negligible. Making 
the conservative assumption that the typical star-formation 
rate remains constant between the observed epoch and the 
present day, the fractional increase in stellar mass from on- 
going star-formation is ~ £lb x (sSFR). Taking the median 
values for the passive galaxy sub-sample implies that the 
fractional increase in stellar mass is only ~ 18% (0.07 dex). 

5.1.1 AGN feedback 

It has been proposed in the literature that mechanical feed- 
back from active galactic nuclei (AGN) could be responsi- 
ble for increasing the size of massive, compact galaxies at 
high redshift (e.g. Fan et al. 2008). In this scenario, AGN 
driven outflows are responsible for evacuating material from 
the central regions of compact galaxies, resulting in an in- 
crease in half-light radius and a decrease in the central 
stellar-velocity dispersion. In fact, because the total mass 



is conserved, the initial and final values of the half-light ra- 
dius and stellar-velocity dispersion are related as follows: 
ojloi = (Rf/R^ 2 . 

Using the dynamical mass measurements for the K- 
Bright galaxies derived here, it is not possible to rule-out 
the AGN feedback mechanism as a viable explanation of 
the observed size evolution. For example, the median stellar 
mass of the K-Bright sub-sample with dynamical-mass mea- 
surements is ~ 10 11 M and the median stellar-velocity 
dispersion is a med = 275 km s _1 . Consequently, if their size 
growth is assumed to be due to AGN feedback, the pre- 
diction would be that SDSS galaxies with stellar masses of 
Mi, ~ 10 11 Mq should have stellar-velocity dispersions of 
a ~ 180 kms -1 . This is consistent with what is actually ob- 
served for early-type galaxies in the SDSS (e.g. van Dokkum, 
Kriek & Franx 2009). 

However, there are several features of the results de- 
rived here which suggest that AGN feedback is not the dom- 
inant mechanism responsible for the size evolution of com- 
pact high-redshift galaxies. The first is our finding that at 
z ~ 1.4 the off-set from the local size-mass relation is not 
a function of stellar population age. The second is that we 
find no indication that the scatter in the galaxy size-mass 
relation is significantly greater at z ~ 1.4 than it is locally. 
Both of these results run counter to the expectations of the 
AGN feedback model as already highlighted by Trujillo et al. 
(2011), based on their study of the size-mass relation of mas- 
sive early-type galaxies at z ~ 1. Moreover, recent numerical 
simulations suggest that the greatest impact of AGN feed- 
back is likely to occur in much younger stellar populations 
than the ^ 1 Gyr old populations seen in massive galaxies 
at z ~ 1.5 (Ragone-Figueroa & Granato 2011). Finally, it 
would also appear unlikely that AGN feedback activity can 
play an important role in compact galaxies with low sSFR, 
given that these galaxies presumably no longer contain a 
large reservoir of gas. 

5.1.2 Major mergers 

Ever since the initial evidence began to emerge that a popu- 
lation of seemingly-old, highly-compact, early- type galaxies 
existed at z ^ 1, the prospect of explaining their subsequent 
evolution via a sequence of major dry (i.e. dissipationless) 
mergers has been widely discussed in the literature. This 
scenario is immediately appealing because the product of a 
major (i.e. mass ratio J?l:3) dry merger is expected to have 
an early-type morphology with a half-light radius which has 
increased proportionally to the increase in stellar mass (e.g. 
Naab et al. 2009). In Fig. 11 the green line illustrates the 
evolutionary path which would be followed by a member of 
the K-Bright sample if its size evolution is driven by major 
mergers. 

Although appealing, there are two obvious problems 
with the major-merger scenario which would appear to ex- 
clude it as the dominant mechanism responsible for the size 
evolution of compact high-redshift galaxies. Firstly, there is 
the problem of the number of major mergers that are ac- 
tually required to reconcile the observed size-mass relations 
at low and high redshift. With reference to Fig. 11, it can 
be seen that to shift the locus occupied by the passive K- 
Bright galaxies at z ~ 1.4 fully onto the local relation size- 
mass relation via major mergers would require a factor of 
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~ 5 — 6 increase in stellar mass. This level of major merging 
is completely inconsistent with N-body simulations, which 
suggest that the typical number of major mergers experi- 
enced in the redshift interval < z < 1.5 by galaxies with 
stellar masses M* 10 11 M is ~ 1 (e.g. Hopkins et al. 
2010). Moreover, this level of stellar-mass growth is totally 
excluded by the constraints provided by the high-mass end 
of the local galaxy stellar mass function. The latest deriva- 
tion of the local galaxy stellar mass function by Baldry et al. 
(2012) suggests that if the passive members of the K-Bright 
sample were to grow in stellar mass by a factor of ~ 5, 
they would overshoot the observed number density of local 
galaxies with stellar masses KU > 1O 11 ' 5 M by a factor of 
~ 25. 

Of course, it could reasonably be argued that the K- 
Bright sample may suffer from significant progenitor bias 
and that the members of the passive sub-sample of the K- 
Bright galaxies are not destined to fully populate the high- 
mass end of the local early-type size-mass relation. If, for 
example, they were only destined to populate the lower half 
of the size-mass relation, the amount of stellar-mass growth 
required would drop to a factor of ~ 2, seemingly making 
the major-merger scenario more plausible. However, if we 
again compare with the local galaxy stellar mass function, 
we find that the effect of progenitor bias on the K-Bright 
sample cannot be severe. Indeed, the results of Baldry et 
al. (2012) show that the number density of local galaxies 
with stellar masses M* > 10 11 M Q is n = (3.4 ± 0.4) x 
10~ 4 Mpc~ 3 . Above this mass threshold, where the K-Bright 
sample should be 100% complete, the number density of K- 
Bright galaxies is n = (3.0 ± 0.5) x 10~ 4 Mpc" 3 (~ 90% of 
the local number density). Even if we restrict the comparison 
to only the passive members of the K-Bright sample, they 
alone can account for ~ 68% of the local number density of 

> 10 11 M galaxies. 



5.1.3 Major plus minor mergers 

Given that there is strong evidence against major mergers 
being able to fully explain the observed size-mass evolution 
between z ~ 1.5 and the present day, it is worth consid- 
ering a toy model in which the additional mass growth is 
dominated by minor mergers, which can increase half-light 
radius in proportion to the square of the fractional increase 
in stellar mass (e.g. Naab et al. 2009). As can be seen from 
the top panel of Fig. 11, this mechanism potentially allows 
rapid growth in half-light radius, without requiring stellar- 
mass growth which violates the constraints provided by the 
local galaxy stellar mass function. 

In the bottom panel of Fig. 11, we illustrate the ex- 
pected evolution of the passive members of the K-Bright 
sample, based on a toy model which incorporates a mix- 
ture of major and minor merging activity. Specifically, the 
toy model assumes that, on average, each passive K-Bright 
galaxy will undergo a single major merger between z ~ 1.4 
and the present day and that the likely mass ratio of this ma- 
jor merger will be 1:3. Additionally, the model also assumes 
that each galaxy undergoes a series of five minor mergers, 
each of which typically has a mass ratio of 1:10. In this 
scenario the stellar mass of the galaxy increases by a fac- 
tor of ~ 2 and the half-light radius increases by a factor of 
~ 3.5. As can be seen in the bottom panel of Fig. 11, this 



~ i — i — i — r~ 



~ i — i — i — r~ 



~ i — i — i — r~ 



n — i — i — i - 




Minor merger _ 
Major merger 



_i i i i_ 



10 



10.5 



11 



11.5 



12 



Log 10 (MyM Q ) 



~ i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r- 




One major merger (1:3) 
Five minor mergers (1:10) 



J I I l_ 



J I I l_ 



10 



10.5 



11 



11.5 



12 



L °gio( M */ M ©) 



Figure 11. The top-panel shows the size-mass relation for the 
passive members of the K-Bright sample, together with the local 
early-type size-mass relation from Shen et al. (2003). The con- 
nected black and green points illustrate the evolutionary path 
followed by galaxies if they undergo minor and major mergers 
respectively. The bottom panel shows the final location of the 
passive K-Bright galaxies if they evolve according to our chosen 
toy model, which consists of a single major merger (mass ratio 
1:3) and five minor (mass ratio 1:10) mergers (see text for a full 
discussion). 



toy scenario is at least capable of reconciling the observed 
size-mass relations for massive passive galaxies at z ~ 1.4 
and the present day, without requiring an unfeasible number 
of major (or minor) mergers. Moreover, in this toy model, 
the descendants of the passive K-Bright galaxies do not vio- 
late the local galaxy stellar mass function because they are 
predicted to contribute ~ 95% of local number density of 
M* > 10 11 M galaxies. 



5.2 Passive disks and the quenching of 
star- formation 

As highlighted in Section 4.2, a substantial fraction of the 
passive members of the K-Bright sample display a disk- 
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like morphology. Indeed, of the forty-one members of the 
K-Bright sample which are classified as passive, accord- 
ing to our adopted sSFR< O.lGyr -1 criterion, nearly half 
(44% ± 12%) have disk-like morphologies (i.e. n < 2.5). We 
note that this is consistent with the recent results of van der 
Wei et al. (2011), who estimated that 65% ± 15% of passive 
galaxies at z ~ 2 with stellar masses M t ^ 6 x 10 10 M© have 
disk-dominated morphologies, based on HST imaging from 
the WFC3/IR Early Release Science (ERS) programme. 

In contrast, the results of a recent study by Bell et 
al. (2012), based on Sersic profile fitting of the CANDELS 
WFC3/IR imaging in the UKIDSS UDS, suggest that at 
z ~ 1.5 there are very few passive galaxies which do not have 
a prominent bulge component. Indeed, Bell et al. (2012) ar- 
gue that a prominent bulge seems to be a necessary, but 
not sufficient, condition for quiescence. Although this would 
appear to contradict the results derived here, it is clearly 
possible that the passive disk-like galaxies identified in this 
study could also harbour a significant stellar bulge, which 
simply cannot be detected using ground-based imaging. 

However, the high fraction of passive disk-like galax- 
ies at z > 1 is confirmed by Bruce et al. (2012), based on 
a detailed disk/bulge decomposition analysis of the CAN- 
DELS WFC3/IR imaging in the UDS. Using a sample of 
~ 200 galaxies with M* ^ 10 11 Mq and photometric red- 
shifts in the range 1.0 < z p hot < 3.0, Bruce et al. estimate 
that 25% ± 6% of passive (sSFR< O.lGyr" 1 ) galaxies are 
disk-dominated (i.e. B/T < 0.5), rising to 40% ± 7% if they 
simply classify disk-dominated objects as having n < 2.5, 
based on a single Sersic fit. It is therefore clear that a sub- 
stantial fraction of massive, passive galaxies at z > 1 have a 
disk-like, or even disk-dominated morphology. 

It is likely that this simple observational fact contains 
crucial information about the evolutionary path followed by 
the progenitors of the massive galaxy population observed 
locally. As discussed in the previous section, the passive 
disks in the K-Bright sample are expected to evolve into 
galaxies with M* ) 2 x IO^Mq, a mass regime which 
is dominated by early-type galaxies in the local Universe. 
Consequently, massive disk galaxies at high-redshift which 
have already been quenched are not a natural prediction 
of evolutionary models in which galaxy-galaxy mergers are 
responsible for quenching star-formation and the morpho- 
logical transformation from late-type to early-type. 

However, the presence of massive, passive disks at 
z ^ 1 would appear to be in good qualitative agree- 
ment with recent hydrodynamical simulations of galaxy for- 
mation. Such simulations suggest that high-redshift star- 
formation is dominated by filamentary inflows of cold gas 
(T < 10 5 K) until a critical redshift of z crit ~ 2 at which 
point, at least for the most massive halos, accretion of hot 
halo gas (T > 10 J K) begins to dominate (e.g. Keres et al. 
2005). Moreover, at z ^ 2 gas is predicted to be shock 
heated to T ~ T vir in halos with masses above a critical 
threshold [Merit — 1O 12 M ), suppressing significant fur- 
ther accretion of gas and efficiently quenching star forma- 
tion (e.g. Dekel & Bimboim 2006). A natural prediction of 
this model is therefore a significant population of massive 
(M* > 1O 1O ' 5 M , M ha i > 10 12 M Q ) quenched disks, which 
are then subsequently transformed into spheroids through 
a combination of both major and minor mergers at lower 
redshifts. 



Finally, we note that this scenario is also in qualitative 
agreement with the phenomenological model of Peng et al. 
(2010), in which star- formation quenching in the highest- 
mass galaxies is independent of environment and propor- 
tional to star-formation rate. This scenario also predicts that 
morphological transformation and star-formation quenching 
are separate physical processes and foresees star-forming 
disks rapidly quenching from the blue cloud onto the red 
sequence. 



6 CONCLUSIONS 

In this paper we have presented the results of an analy- 
sis of the relationship between galaxy size, stellar mass and 
specific star-formation rate in a mass complete (M* 6 x 
10 n Mo) sample of spectroscopically confirmed 1.3 ^ z ^ 
1.5 galaxies in the UKIDSS UDS field. Using stellar-mass 
measurements based on spectro-photometric SED fitting 
and galaxy size measurements based on two-dimensional 
modelling of the available UDS K— band imaging, we have 
investigated the evolution of massive galaxies on the size- 
mass plane as a function of specific star-formation rate, 
galaxy morphology and stellar population age. In addition, 
we have derived an alternative constraint on the average size 
evolution of massive galaxies since z ~ 1.4 using dynami- 
cal mass measurements for a small sub-sample of galaxies 
with reliable stellar-velocity dispersion measurements. Fi- 
nally, we have investigated evolutionary scenarios which can 
plausibly explain the observed size evolution without vio- 
lating the constraints provided by both the latest N-body 
dark-matter simulations and current determinations of the 
local galaxy stellar mass function. The main conclusions of 
this study can be summarized as follows: 

(i) The location of massive galaxies (M > 6 x 10 10 M Q ) 
on the size-mass plane at z ~ 1.4 is more closely correlated 
with sSFR than galaxy morphology. 

(ii) Massive galaxies at z ~ 1.4 which are passive (sSFR 
< 0.1 Gyr" 1 ) follow a tight size-mass relationship, which is 
identical in slope (and scatter) to that displayed by local 
early- type galaxies in the SDSS. However, at a given stel- 
lar mass, passive galaxies at z ~ 1.4 have half-light radii 
a factor of f g = 2.43 ± 0.20 smaller than their low-redshift 
counterparts. 

(iii) In contrast, massive star- forming galaxies at z ~ 1.4 
(sSFR > 0.1 Gyr -1 ) lie closer to the size-mass relation of 
local late-type galaxies, but are still a factor of f g = 1.61 ± 
0.17 smaller than their low-redshift counterparts at a given 
stellar mass. 

(iv) If the K-Bright sample is split into early and late- 
type morphologies using a Sersic index threshold of n = 2.5, 
the corresponding growth factors are found to be: f g = 
2.37±0.29 and f g — 2.15±0.15 for early and late-type galax- 
ies respectively. However, the early-type galaxy size-mass re- 
lation at z ~ 1.4 displays a scatter which is noticeably larger 
than that associated with the passive galaxy size-mass re- 
lation. Moreover, the correspondence between passivity and 
early-type morphology (and active star-formation and late- 
type morphology) in the K-Bright sample is found to be 
weak. 

(v) Using a small sub-sample of the K-bright galaxies 
with reliable stellar-velocity dispersion measurements, an al- 
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ternative constraint on the average size evolution between 
z ~ 1.4 and the present day is obtained. The derived value 
of f g = 2.33 ± 0.32 is entirely consistent with our primary 
determination based on stellar-mass measurements. 

(vi) We demonstrate that a toy model which combines a 
single major merger with a sequence of five minor merg- 
ers can reproduce the observed size evolution of massive 
galaxies, without violating constraints imposed by the local 
galaxy stellar mass function and the predictions of N-body 
simulations. 

(vii) The significant fraction of passive galaxies with disk- 
like morphologies in the K-Bright sample provides additional 
evidence that separate physical processes are responsible for 
the quenching of star-formation and the subsequent mor- 
phological transformation in massive galaxies at z ~ 1.4. 
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APPENDIX A: THE K-BRIGHT SAMPLE 

In Table Al we provide the basic observational properties 
and important derived quantities for the eighty-one objects 
in the K-bright sample. 
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Table Al. Basic measured and derived properties of the K-bright sample. The first three columns list the catalogue ID number and 
coordinates of each member of the sample. Columns four and five list the spectroscopic redshift (typical uncertainty <5z = 0.001) and 
total K— band AB magnitude (uncertainty ^ 5%) for each object. Columns six and seven list the derived half-light radii (circularized) 
and stellar masses along with their respective uncertainties. Columns 8-10 list the Scrsic index (typical uncertainty 30%), specific star- 
formation rate and stellar population age (typical uncertainty 0.2 dex) for each object respectively. The typical uncertainty on the derived 
sSFR values is dominated by the scatter associated with the adopted star-formation calibrations and is a minimum of a factor of two. 
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0.54 
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The galaxy size-mass relation at z ~ 1.4 

Table Al. Continued. 
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